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Abstract 

We show that it is possible to obtain self-consistent and physically accept- 
able relativistic classical equations of motion for a point-like spin-half particle 
possessing an electric charge and a magnetic dipole moment, directly from a 
manifestly covariant Lagrangian, if the classical degrees of freedom are ap- 
propriately chosen. It is shown that the equations obtained encompass the 
well-tested Lorentz force and Thomas-Bargmann-Michel-Telegdi spin equa- 
tions, as well as providing a definite specification of the classical magnetic 
dipole force, whose exact form has been the subject of recent debate. Ra- 
diation reaction — the force and torque on an accelerated particle due to its 
self-interaction — is neglected at this stage. 



I. INTRODUCTION 

The "classical limit" of mechanics has always played an important role in practical 
physics. While it may be regarded as preferable to describe the behaviour of a physical 
system by solving the problem exactly within quantum mechanics, in practice one need 
not always do so. If one is only interested in expectation values of operators, then one can 
make good use of Ehrenfest's theorem — which holds good under quite general circumstances 
(see, e.g., §31) — and instead solve the equations of classical Hamiltonian or Lagrangian 
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dynamics. Thus, for example, the Lorentz force law for charged particles, and the Thomas- 
Bargmann-Michel-Telegdi spin precession equation for particles with spin, are used to 
advantage every day — implicitly or explicitly — in a wide variety of practical situations. 

These two examples, however, owe their widespread acceptance largely to an important 
quality they possess: they can be derived simply from a knowledge of the external elec- 
tromagnetic multipolar fields of their point-like sources, regardless of the detailed structure 
of these sources On the other hand, the force on a point-like particle with spin is not 
uniquely defined by the external properties of its dipolar electromagnetic fields; this fact, 
well noted by Thomas in 1927 |@], and more recently by Hrasko 0], is not yet universally 
appreciated. To obtain the "correct" force on (and, hence, complete equations of motion for) 
a particle with a magnetic moment, one must therefore make some assumptions, explicit or 
otherwise, about the nature of the object creating the dipolar field; it is this freedom that 
has contributed to the ongoing controversy on the subject 0-0. 

It would be impossible to review here all of the various assumptions underlying previous 
attempts at a complete classical description of particles with spin; instead, we shall simply 
state the premises and results of a theoretical analysis that we have carried out. In section 
we outline our considerations in choosing a suitable relativistic classical Lagrangian for a 



particle with spin. We then, in section |T|, outline our derivation of the equations of motion 
from this Lagrangian, and present the results in what we believe are the most transparent 
forms. As will be seen, the technique we shall use leads to the omission of all aspects 
of radiation reaction — which is an omission that should be repaired — but nevertheless the 
results obtained do allow extensive contact with existing knowledge about the classical 
behaviour of particles with spin. 

II. AN APPROPRIATE LAGRANGIAN 

As already mentioned, some non-trivial assumptions must be made for one to obtain 
the force equation of motion for even a point-like particle endowed with a magnetic dipole 
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moment. Our approach is to use a Lagrangian approach from the outset, and hence ben- 
efit from the abovementioned guarantee that Ehrenfest's theorem provides. This route is, 



however, still a potential minefield; as recently emphasised by Barut and Unal |jTO[, the 
"spinning top" (or "current loop") phase space degrees of freedom are not appropriate to 
quantum spin-half particles. We shall return to this question shortly. 

Our first consideration is to outline precisely what type of system we wish to address. In 
this specification we shall be exacting: we shall only be considering the intrinsic magnetic 
dipole moment of a spin-half, pointlike, structureless particle — together with any electric 
charge interaction it possesses, of course — in externally applied (but, in general, time- and 
space-varying) electromagnetic fields. It is only with such a tight restriction of our focus 
that we can proceed confidently in the classical realm at all, as will become apparent. For 
definiteness, we consider the electron a good approximation to the type of particle we are 
considering; however, we shall not claim any predictive control over contributions from the 
anomalous magnetic moment of the real electron, arising as it does from QED and other 
processes that are, in reality, additional degrees of freedom over those properly described by 
the single-particle Dirac equation. We shall not enter the debate as to whether a particle 
might "intrinsically" have a g-factor differing from 2 in the single-particle Dirac equation; if 
so, then such "anomalous" moments would be included in our considerations. In any case, 
since the pure Dirac moment of the electron dominates its anomalous moment numerically, 
our considerations will, at least, provide useful practical results for real electrons to leading 
order in the fine-structure constant. 

Our main task in this section is to explain the particular Lagrangian that we have chosen 
to represent the magnetic interaction of such a particle. There are two distinct parts to this 
decision: firstly, the choice of the value that the Lagrangian should take; and, secondly, 
a determination of the appropriate classical degrees of freedom that are to be used in the 
Euler-Lagrange equations — and, hence, the functional form of the Lagrangian. The former 
choice is laid out for us: the Dirac equation gives us 
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Lint = -|^^°^F./3, (1) 

where q is the charge of the particle (= — e for an electron) and g is its gyromagnetic ratio 
(= 2 for the simple minimally-coupled Dirac equation), and our units follow SI conventions 
with the exception that we set £o = A*o = c = 1. However, to entertain the possibility 
of the single-particle Dirac equation fully describing a neutral particle with an anomalous 
moment, we shall replace the g factor by the corresponding magnetic moment fi = gqh/im, 
thereby separating the interaction of electric charge from that of magnetic moment, at least 
formally. 

Our second task is infinitely more perilous: choosing an appropriate set of generalised 
coordinates to represent a spin-half particle. Firstly, we observe that the position of a 
particle is usually assumed to be an appropriate quantity in the classical limit; we shall also 
make this assumption. Thus, we take the (expectation value of the) four-position z°'{t) to 
be four of the independent degrees of freedom. (One can show that, in a general dynamical 
framework allowing particles with time-dependent masses, all four of these coordinates are 
indeed independent.) If we take the proper time for the particle, r, to be the generalised 
time of the classical Lagrangian framework, then the generalised velocities corresponding to 
the z"" are given by 

[/" = ^ = z". (2) 
or 

Having defined the four- velocity thus, we can then simply define the centre of energy frame 
for the particle as that instantaneous Lorentz frame in which = (1, 0, 0, 0). 

We now turn to the magnetic moment of the particle. For a spin-half particle, the 
magnetic moment is parallel to the spin. In the centre of energy frame, we expect the 
particle's spin to be describable by a two-component spinor, of fixed magnitude s = \Ti. 
The expectation value of this spinor is equivalent to a fixed-length three-vector , in this 
frame, whose direction describes simultaneously the expectation value of polarisation along 
any arbitrary axis, as well as the phase difference between the "up" and "down" components 
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along such an axis, but which throws away the overall phase of the wavefunction (as does 
any reduction to the classical limit). It is this fixed-length three- vector that we shall assume 
to represent the magnetic moment of the particle. 

Consider now the dynamical nature of spin angular momentum in the classical limit. 
It would seem, from the previous paragraph, that the spin should also be considered as 
a fixed-length three-vector in the centre of energy frame of the particle. However, this 
would preclude writing down the usual classical kinetic Lagrangian of rotational motion in 
terms of generalised velocities {i.e. the time derivatives of the Euler angles), namely, in the 



non-relativistic limit, ITT 



L = \uj-I-uj, (3) 

and thereby preclude obtaining the torque on the particle by means of the Euler-Lagrange 
equations for the rotational degrees of freedom. Our approach, therefore, shall be the fol- 
lowing: we shall consider the spin of the particle to be given, in the non-relativistic limit, 
by its usual form 



s = I ■ UJ 



with the "kinetic" term (§) retained in the Lagrangian in this same limit. This allows the 
magnitude of the spin, \s\ = s to be an a priori dynamical variable. (We know, of course, that 
s should, ultimately, be a constant of the motion {i.e. ^h) if we are to accept the equations 
of motion as applicable to spin-half particles; we shall discuss this necessity in greater detail 
shortly.) We shall, on the other hand, still consider the magnetic moment of the particle to 
be formally a function of the generalised coordinates themselves {i.e. the Euler angles) and 
not the velocities. This distinction between the spin and magnetic moment of an electron is 
an unfamiliar concept within quantum mechanics; indeed, in that situation, no distinction 
need be made. However, this distinction is vital in the classical framework — at least, if one 
wants to derive the equations of motion from a Lagrangian — as the spin angular momentum 
and magnetic dipole moment of an arbitrary classical object are not, in general, related in 



any special way. Our ansatz that the magnetic moment is dependent on the Euler angles 
alone cannot be justified a priori; the only justification will be that it produces equations 
of motion that satisify the strict requirements of consistency for spin-half particles — that 
have not, to our knowledge, been satisfied by previous classical Lagrangian approaches. We 
further note, at this stage, that the ultimate parallelism of spin and magnetic moment does 
not bar us from considering them to have different functional dependencies: the fact that 
they are parallel may be introduced, in this classical context, as simply a "constitutive" 
relation, which can only be fully justified upon investigation of the quantum mechanical 
analysis. 

We now turn to synthesising the information that we have laid out above in a self- 
consistent, relativistic, Lagrangian framework. Firstly, it is necessary to generalise the 
non-relativistic concept of the spin angular momentum s to the relativistic domain. This 



procedure is carried out in many textbooks on electrodynamics (see, e.g., [0, pp. 556- 
560); we shall not add anything new. One simply defines a four- vector 5" such that, in 
the rest frame of the particle, it has vanishing zero-component S^, and its three-vector part 
S is equal to the non-relativistic spin s. Since, by definition, the three-vector part of the 
four-velocity, U, vanishes in this frame, the identity 

S"Ua = (4) 

shows that only three of the 5" are independent, as one would expect from the non- 
relativistic case. 

In a completely analogous way, one can generalise the non-relativistic angular velocity 
vector (jj and the moment of inertia tensor / to their relativistic counterparts cj" and I"^. 
The non-relativistic rotational "kinetic" Lagrangian (j^) can then be written relativistically 

as 



^rot 



^CU'^IafSUJ^ (5) 



One might wonder, at this point, how a spin-half particle can have a "kinetic" term of 
rotation, when it is well known that there is no classical "rotating model" that represents 
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spin angular momentum. The answer is as follows: the functional form of the spin kinetic 
Lagrangian must be retained in the classical limit, regardless of whether or not it corresponds 
to any particular "model" that one might dream up. In fact, it will be found that the 
quantities a;° and will completely disappear from the final equations of motion; the only 
remaining quantity will be the physically observable quantity S"'. Thus, by what might 
arguably be considered a sleight of hand, we can analyse the system in question within 
Lagrangian mechanics, without recourse to the particular Newtonian models of refs , 
nor the original heuristic (albeit brilliant) arguments of Thomas and Bargmann, Michel 



Before we can write down the final expression for the Lagrangian we shall use, we must 
first "massage" the magnetic interaction Lagrangian (||) into a more suitable form. Our 
major task is to the interpret the spin tensor, cr"^, that the Dirac equation introduces. This 
is not a trivial task: translating between the spin vector that we have already defined, S*", 
and a spin tensor, requires the use of both the alternating tensor £°'f^t^'^^ and another four- 
vector. Often, in quantum mechanics, one uses the canonical momentum vector, p"", for this 
purpose; the resultant "spin" vector is known as the Pauli-Lubanski vector. 



This quantity is, indeed, very useful in many situations. However, we are here considering 
Lagrangian mechanics; therefore, we should expect that mechanical momenta (or, in other 
words, generalised velocities) should be employed; canonical momenta belong to Hamiltonian 
dynamics. We therefore follow Jackson ( |T^, p. 556) in using the four-velocity to define the 
transformation between the spin tensor and the spin vector, namely 



and Telegdi [^]. 



(6) 



'a 



(7) 



It is straightforward to verify that the reverse transformation of (|^ is given by 



S' 



'al3 ^a/Bfiujj- ^ 



(8) 
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Insering this into (0) and simplifying, we finally obtain our desired magnetic interaction 
Lagrangian, 

Ur.t = /i"F„^f/^ (9) 

where F^p = \sai3^'^F^i, is the dual electromagnetic field strength tensor, and the magnetic 
moment four-vector, /x" is, as noted previously, considered to be a four-vector "embedded" 
in the instrinsic rotational coordinates of the particle, and hence is functionally dependent 
on the Euler angles, but not their derivatives. 

Our complete Lagrangian is then assembled simply from the kinetic rotational term (^, 
the magnetic interaction term (^), and the standard translational kinetic and electric charge 
interaction terms: 

L = \mU''U^ + iu;"J„^ ujf^ + g[/°A„ + fi'^F^pU'^. (10) 

It is this Lagrangian, and most particularly its functional form, that forms the basis of the 
following section. 



III. DERIVATION 

We now turn to the question of deriving, from the Euler-Lagrange equations, the equa- 
tions of motion for the particle under study. (In the following, our language shall describe 
classical quantities in the same way that Newton would have done, but in reality we are 
referring to the expectation values of the corresponding quantum mechanical operators for 
the single particle in question.) The generalised coordinates for the particle, following the 
discussion of the previous section, are taken to be the four translational degrees of freedom 
2;", together with the three Euler angles describing the intrinsic "orientation" of the particle. 

The mathematical manipulations necessary to obtain the seven Euler-Lagrange equa- 



tions are in principle no different to those in everyday classical mechanical problems [ITT 



There is one subtlety, however, that enters into one's consideration of the spin vector S'^ 
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and the magnetic moment vector /x". By their very definition, these vectors are "tied" to 
the particle's four-velocity, in the sense that identities such as (^) always holds true. One 
must therefore be careful when defining their proper-time derivative: since they are, in ef- 
fect, defined with respect to an accelerated frame of reference, there is a difference between 
taking the time-derivatives of the components of the vector, and the components of the 
time-derivatives of the vector, as General Relativity teaches us. In fact, the philosophical 
framework of General Relativity tells us that it is the latter that is the generally invariant, 
"covariant" derivative which should be used in the relativistic Euler-Lagrange equations; 
the former is simply the "partial" derivative. However, it is straightforward to verify that 
they can be related via 

d ^(ja^^a (U^C^) ^ (11) 

where C" is any space-like vector that is orthogonal to the particle's four-velocity (such 
as 5" and fi"), the left hand side denotes the "covariant" derivative, and the "partial" 
derivative, of such a vector. (The Thomas precession pr5|,3 is, in fact, just another way 



of expressing this "General Relativity" effect — namely, the non-vanishing commutator of 
Lorentz boosts.) 

It is now relatively straightforward to apply the Euler-Lagrange equations to the La- 
grangian (0). The three Euler angle degrees of freedom lead, just as in the non- relativistic 
1TT| , to the torque equation of motion for the particle. With some algebra and simpli- 



case 



fication, they (together with the identity (^ and the relation ([Tl|)) lead to the four- vector 
equation of motion 

+ U^if^Sp = F^pfi^ + U^fi'F.pU". (12) 

It should come as no surprise that this equation is precisely that obtained by Bargmann, 
Michel and Telegdi (equation (6) of ref. 0), before they simply substituted in the Lorentz 
force law for f/g {i.e. equation (7) of [Q). As was noted in section |I|, this equation — (|T^) — is 
true in general for any point-like object generating an external magnetic dipole field 0]. 



We now turn our attention to the four translational degrees of freedom of the parti- 
cle. The Euler-Lagrange equations for these coordinates, for the same Lagrangian (|10|) 
as used above to generate the Bargmann-Michel-Telegdi equation, immediately yield the 
four-equation of motion 

The first term on the right side is, of course, simply the Lorentz force. The second and third 
terms, on the other hand, are infinitely more interesting. The second term represents the 
gradient forces on the magnetic dipole moment, i.e., in the non-relativistic limit, (/x • V) B. 
The third term represents a force of the type — /x x E in the non-relativistic limit. The 
reason that they are so interesting is that many authors have argued that precisely these 
terms should constitute the non-relativistic limit of the magnetic dipole force (see, e.g., 
and references therein; also Here, we have obtained these terms from a relativistic 

Lagrangian directly, without need for assumptions other than those outlined in the previous 
section. 

The last term in (|T3]) is, as a matter of principle, much deeper, but in practice completely 
ignorable. It constitutes a contact force between the particle and the electromagnetic current 
that is generating the "external" field. In practice, such an interaction is usually ignored; 
however, upon further investigation, it is recognised that it is precisely this contact force 
that allows the expression (|13]) to otherwise resemble so closely the force on a "monopole- 
constructed" dipole model 0,0, while on the other hand possessing an interaction energy 
equivalent to a "current loop" model @J^-§| as required for agreement with atomic hyperfine 
levels |]T^. It is in this way that the "Lagrangian-based" model of this paper seems to select 
those desirable properties of both the "monopole-constructed" and "current-loop" models, 
while being classically equivalent to neither. 

Now that we have dealt with the dynamics of the Euler-Lagrange formalism, and have 
obtained allegedly appropriate equations of motion, we can now specify the exact "constitu- 
tive relation" between the magnetic moment and the spin of the particle, for the particular 
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case we wish to study: a spin-half particle. From quantum mechanics, we know that the 
these two quantities are always parallel, namely, 

= (14) 

where ( = fi/s is some constant, a property of the particle in question, which, for particle 
of charge q, is commonly written in terms of the g factor as C = gq/2m. Before we can use 
(|l^), however, we first note that equation (|13|) is itself in a somewhat awkward form: the 



right hand side involves both /i" and ?7". However, the parallelism condition ([14|) allows us 



to uncouple equations (]1^) and (|13]) simply by substituting ([1^) into (|T3D, yielding 



^ (m[/„) = qF^pU^ + U^fi'd^F^p + ( {F^pF^'^fi'' + F^pU" {fi'^F^Mn} , (15) 
where we have, for practical simplicity, dropped the contact term in ([T3|). 



There are now several important things we can do with equations (0) and (|T^). Most 
importantly, we examine how the mass of the particle changes with time; if it is not constant, 
then our equations of motion cannot possibly apply to any particle, such as an electron or 
muon, that has a constant rest mass. (It should be noted that the classical relativistic 
formalism we have employed has made no assumptions as to the constancy of the rest mass 
m with time: in general, the mass may change as the system in question gains energy from or 
loses energy to the external field.) It is straightforward to verify that the general proper-time 
rate of change of the mass of a particle, rh, can be expressed as 

m = U''^{mUa). (16) 

(Jj I 

Using the the identity (^), and the fact that FapF^^ is proportional to the metric, gav^ it 
follows that ( p!5D yields, in (|16|), m = 0. In other words, the set of equations (|T2D and (|I5|) 



rigorously maintain the constancy of the mass of the particle. This property is far from 
trivial; for example, in a classic textbook ( |T5[, p. 74), a time- varying "effective mass" of 
the electron was introduced, due to the fact that the equations of motion derived therein 
allowed time-changing rest masses. It is of utmost importance that this difficulty — present 
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in most previous attempts at consistent equations of motion for a dipole — is overcome in 
the Lagrangian treatment of this paper. 

We now verify that the the parallehsm condition, (0), is consistent with the equations of 
motion (p!^ ) and (|15|). The reason for this concern is that we have assumed that the magnetic 
moment is a const ant- magnitude vector, whereas the spin itself is a dynamical quantity that 
may, in general, change its magnitude. We must verify that, if (|l^ is assumed to hold true 
at one particular proper time r, the equation of motion ([T^) does not change the magnitude 



of the spin. On diffentiation of the identity = —S'^Sa, one finds s = —^S°'Sa', use of ( p!4D 
in (|12D then shows that (|l^) does, in fact, yield s = 0. Thus, the magnitude of the spin is a 
rigorous constant of the motion, just as is the mass. 

Our work is now essentially complete. However, to make practical use of the equations 
(|l^) and (p!5D, it is appropriate to both present them in a more computationally- friendly 
form, and to highlight clearly where they add to existing knowledge of spin-half particles. 
These two tasks can essentially be carried out simultaneously. We shall transform ([T^ ) and 
(p^) into the standard "3 + 1" form, as, for example, presented in Jackson's textbook ( ||12|| , 
pp. 556-560). This involves using the three- velocity of the particle in some particular fixed 
"lab" frame, v (Jackson uses the symbol (3), as well as the three-spin s of the particle as 
seen in its rest frame, but referred to the (non-rotating) coordinates of the "lab" frame. 



The procedure used to effect this transformation is described in detail in |]12[; the algebra is 
tedious, but straightforward. The results for equations (|15D and (|l^) are 

at 7m 7m 7^m 



and 



— = sx flnew, (18) 



where 



[ 7 m ) [ 7-l-lmJ 7 + 1 I m) 

H if/ H V X B [v -Ejv ^ -V X s 

7 m 7 + lm 7 + lm m(7 + 1) 
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and we have defined the convenient quantities 



s 

ot 7 + 1 

s = s — [v ■s)v, 

7 + 1 

E' = E + vx B — (v -E) V, 

7 + 1 

E" = E + vx B -{v ■E)v, 
B" = B -vx E-{v ■B)v, 
and (?eff = C ■ V) + 7C -v) & - 7C' {s -E') , 

and, in all expressions, the partial derivatives act only on the external field quantities E 
and B. 

For ease of comparison with the equations of motion in current usage, we present the 



Lorentz force law in the same form as (|T^, 



= ^E", (19) 

at 7m 

and, likewise, the precession frequency vector for the Thomas spin equation ( |12|, p. 559): 

[ 7mJ [ 7 + lmJ 7 + II m) 

It can be seen that, as advertised, the Lorentz and Thomas equations are contained com- 
pletely in the new equations. However, several new features are present in both the new 
force equation, ([T^), and the new precession frequency vector, flnew Most obviously, the 
magnetic dipole force is now included in (|r^), albeit somewhat obscured by the multitude of 
"convenient quantities" introduced for typographical sanity. A recognition of this expression 
may again be obtained by taking the non-relativistic limit (first order in v, ignoring Thomas 
precession); ( [17|) then returns us to 

^(mv) = q (E + V X B) + {fi - V) (B - V X E) ~ /l X E, 
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which is, as noted earher, the now generally-accepted 0,1^-3 dipole force expression. There 
are, of course, numerous new subtleties of ( ]T7| ) that arise from relativistic kinematics; we 
shall, however, defer a more exhaustive investigation of them to another place. 

Turning, now, to the spin precession equation (|1^), what may come as a surprise to some 
is that there are differences between floid and finew The reason, however, is simply found: 
the Thomas precession effects in the spin precession equation depend on the acceleration 
of the particle; clearly, if we have obtained a more accurate force law, then these changes 
will necessarily feed through to the spin equation as well. It is interesting to note that 
the approximate nature of the Thomas expression floid was noted explicitly by Thomas 
himself 0, and by Bargmann, Michel and Telegdi in their re-derivation ||^, but that, in 
the intervening decades, this approximate nature has been lost on many practitioners, who 
mistakenly believe fioia to be an exact expression. 

At this point, it is worthwhile commenting on the absence of any presence of radiation 
reaction forces in the equations that have been derived above. Their omission can, in fact, 
be traced back directly to the (implicit) assumption that the electromagnetic potentials and 
fields in the Lagrangian (|T^) are the externally generated ones only. In the case of the 
electric charge interaction, it was recognised by Lorentz []TB[, and later clearly explained by 



Heitler [|17|, that this is an incorrect assumption. For full consistency, one must include the 
potentials and fields of the particle itself in the Lagrangian, even though they appear, at 
first sight, to be hopelessly divergent quantities. If one proceeds with extreme care, one 
can show that the effects on the equations of motion of this "self-interaction" are, in fact, 
twofold: firstly, the addition of the (finite) radiation reaction terms of the Lorentz-Dirac 
equation; and, secondly, the dynamical effects explaining the electromagnetic self-energy 
contributions to the rest mass of the system, for whatever arbitrary charge distribution is 
assumed. The beauty of this procedure is that it reveals the all-encompassing nature of the 
Lagrangian; no further input is required to obtain the equations of motion. 

The success of the above procedure clearly indicates that a similar procedure should 



be undertaken for the additional magnetic moment interaction Lagrangian present in (|10|). 
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Barut and Unal [|T0| have considered the dipole radiation reaction question from the point 
of view of a semi-classical Zitterbewegung model, but to our knowledge an exact treatment 
of this problem in terms of the classical spin vector s has not been performed. This problem 
is one that we are currently investigating. However, it is vastly more complicated than the 
electric charge case, by virtue of the inclusion of rotatonal degrees of freedom for the particle. 
In addition, one already knows in advance that the electric charge and magnetic moment 
interact in their respective radiative terms, as is evidenced by the Sokolov-Ternov effect 
|T8| (the polarisation of electrons due to their emitted synchrotron radiation), spectacularly 
confirmed in the polarisation experiments at LEP in recent years [|19|,^. Any prospective 
solution of the combined radiation reaction equations for a charged, spinning particle must 
therefore, at the very least, reproduce this important result. 

IV. CONCLUSIONS 

It has been shown, in this paper, that it is possible to construct a fully consistent, 
comprehensive, relativistic classical Lagrangian framework for analysing the motion of spin- 
ning particles possessing both electric charge and magnetic dipole moments. The results 
obtained here are not revolutionary. They encompass the well-known Lorentz force and 
Thomas-Bargmann-Michel-Telegdi equations. They provide a rigorous foundation for the 
magnetic dipole force law currently believed to be the most appropriate for such particles. 
They further integrate, seamlessly, this dipole force with the Lorentz force and Thomas- 
Bargmann-Michel-Telegdi equations, in a fully relativistic way. 

It should be noted that the results of this paper agree with the lowest-order terms 



obtained in the analysis of Anandan, based on the Dirac equation ||2T| — and, of particular 
note, the "Anandan force" proportional to E x [fi x B), which is, of course, simply the 
zeroth-order term in the force term —fi x E above, when it realised that /x is, via the Thomas- 
BMT equation, proportional to /x x S to zeroth order. The current work, however, includes 
terms to all orders in the particle's velocity, not just the lowest-order limit of Anandan. It 
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should also be noted that the criticisms of Casella and Werner of Anandan's analysis 



21 1 are erroneous, being based on an obvious omission of all "spin- flip" terms from their 



quantum mechanical equations of motion. 

The framework outlined in this paper may now be used as a platform for full inclusion of 
radiation reaction — not just for the electric charge, but also for the magnetic moment, and 
their mutual interactions — in the classical limit. If the enormous assistance provided by the 
existing classical radiation reaction theory, both in terrestrial and astrophysical applications, 
is any guide, then one can only speculate what additional richness of physical phenomena 
will be made sensible with this addition to our analytical resources. 
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